% 输入参数
Re = 185.1;         % 定子电阻
Xm = 9.9554;        % 磁化电抗

R1 = 0.2268;        % 一次绕组电阻
R2pr = 0.12528;     % 二次绕组电阻（参考值）
X1 = 0.51252;       % 一次绕组电抗
X2pr = 0.76878;     % 二次绕组电抗（参考值）

Rc = 185.1;         % 异步电动机定子电阻
Xm = 9.8554;        % 异步电动机磁化电抗

V1 = 230/sqrt(3);   % 一次电压（线电压）
f = 60;             % 频率
p = 6;              % 极对数

Pfw = 224;          % 定子铜损耗
n = 2.8;            % 定子效率

npts = 200;         % 绘图点数
s = linspace(0.00001, 1, npts);    % 定义s值
s = fliplr(s);      % 反转s值数组，从大到小排序

% 初始化变量
I1 = zeros(1, npts);    % 输入电流
TTd = I1;               % 发展转矩
PF = I1;                % 输入功率因数
Ps = I1;                % 输出功率
eff = I1;               % 效率
nm = I1;                % 机械转速

ws = 2/p * 2*pi*f;      % 同步角速度
ns = 120*f/p;           % 同步转速

R2pr0 = R2pr;           % 保存初始R2pr的值
X2pr0 = X2pr;           % 保存初始X2pr的值
smax = R2pr/sqrt(R1^2 + (X1 + X2pr)^2);    % 计算最大s值

% 循环计算感应电机性能
for i = 1:npts
    Z2 = R2pr/s(i) + j*X2pr;
    Zm = j*Rc*Xm/(Rc + j*Xm);
    Zin = R1 + j*X1 + Z2*Zm/(Z2 + Zm);
    I11 = V1/Zin;
    I1(i) = abs(I11);
    PF(i) = cos(angle(I11));
    I2pr = abs(Zm/(R2pr/s(i) + j*X2pr + Zm)*I11);
    TPin = 3*V1*I1(1)*PF(i);
    TTd(i) = 3*I2pr^2*R2pr/s(i)/ws;
    nm(i) = (1 - s(i))*ns;
    PPo = TTd(i)*(1 - s(i))*ws - Pfw*(nm(i)/ns)^n;
    if PPo < 0
        break;
    else
        Ps(i) = PPo;
        eff(i) = 100*Ps(i)/TPin;
    end
end

% 绘制图形
subplot(2,1,1),
plot(nm, TTd);
grid;
title('Developed torque');
xlabel('Speed, rpm');
ylabel('Torque, N-m');

subplot(2,1,2),
plot(nm, Ps/746);
grid;
xlabel('Speed, rpm');
title('Output power');
ylabel('Output power, hp');

figure(2);
subplot(2,1,1);
plot(nm, I1);
grid;
title('Input current');
xlabel('Speed, rpm');
ylabel('Current, A');

subplot(2,1,2);
plot(nm, eff);
grid;
title('Efficiency');
xlabel('Speed, rpm');
ylabel('Efficiency, %');

figure(3);
subplot(2,1,1);
plot(nm, PF);
grid;
title('Input power factor');
xlabel('Speed, rpm');
ylabel('Power factor');